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Abstract 

We investigate numerically using the bond-fluctuation model the adsorption of a random 
AB-copolymer at the interface between two solvents. From our results we infer several scal- 
ing relations: the radius of gyration of the copolymer in the direction perpendicular to the 
interface (Rgz) scales with the interfacial selectivity strength, as Rgz — N'^ f{-\/Nx) where 
V is the usual Flory exponent and is the copolymer's length; furthermore the monomer 
density at the interface scales as x^" for small x- We also determine numerically the monomer 
densities in the two solvents and discuss their dependence on the distance from the interface. 
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Copolymers at interfaces are very important in technical applications. For example, the interface 
between two immiscible polymer melts can be mechanically reinforced by dissolving copolymers 
into the system [Q, In other applications the surface tension between two immiscible solvents, 
e.g., oil and water, can be reduced by adding block copolymers consisting of a hydrophilic and a 
hydrophobic part. Thus the adsorption of copolymers at interfaces has received particular atten- 
tion [3-7]. For diblock copolymers, the difference in the solubilities of the monomers favors the 
localization of the copolymer at the interface, with each block in its favorable solvent. However, 
for random copolymers frustrated situations may arise since the chain's connectivity forces some 
monomers to stay in their unfavourable solvent. Using a Hartree-type approach Garel et al |3j 
have studied the localization transition of an ideal random chain at an interface. Yeung, Balazs 
and Jasnow Q have addressed the question of correlations in the A- and B-distributions, a point 
investigated recently through computer simulations by Balazs et al j^] in the context of copolymer 
brushes. In this paper we study numerically under excluded volume conditions a single random 
copolymer at an interface and we pay particular attention to scaling . 

We use the bond-fluctuation model (BFM) |l^ to perform Monte Carlo (MC) simulations. In 
the BFM the polymers obey the excluded- volume requirements, and the motion occurs such that 
the bonds do not cross each other, see Refs. |l^ for details. Here we take a cubic box of size 
L X L X H , with periodic boundary conditions in the x- and y-directions and two impenetrable 
surfaces at z = and z = H. We study the behavior of a single copolymer consisting of N ran- 
domly chosen monomers of A- and of B-type. We assume a symmetrical situation: the interaction 
parameter of the monomers is x^bT when immersed in their unfavorable solvent and zero other- 
wise. The solvent below the interface (z < H/2) favors A-type monomers, and the solvent above 
the interface (z > H/2 + 1) B-type monomers. Note that the interface is thus at z = {H + l)/2. 
In each Monte Carlo step the chain moves by position changes of the monomers, which attempt 
nearest neighbor steps on the underlying cubic lattice. A move is taken into consideration only if 
it satisfies the requirements of self-avoidance and of non-crossing of bonds. Furthermore, energet- 
ically unfavorable moves are statistically permitted according to the usual Boltzmann factor. 



We obtain results for copolymers of lengths N = 16, 32, 64, and 128 using systems with sizes 



L = 50 and H = 100. An initial configuration is generated starting with the first monomer near 
the interface and then randomly adding the subsequent monomers such that self-avoidance and 
non-crossing of bonds are obeyed. The energetic aspects of the interaction with the solvents are 
then taken care of by the usual Boltzmann factor; the monomer-monomer interaction is only ac- 
counted for through the excluded volume aspects. We let the chain move for a long time according 
to the MC prescriptions, such that the chain relaxes to equilibrium. The averaged quantities which 
we will show below are then obtained from such equilibrium configurations. We found numerically 
the relaxation time (determined using the autocorrelation function of the radius of gyration Rg 
and of its z~component Rg^ Q) to be around 50,000 Monte Carlo steps (MCS), where a MCS 
consists of N move attempts; we thus view the copolymer as having reached equilibrium after 
200,000 MCS. Averages are then calculated from the configurations obtained in the subsequent 
200,000 MCS. For each copolymer length we realized 100 independent runs. 

Figure 1 shows the normalized probability to find monomers of A-type and of B-type at the height 
z. In this example the copolymer length is = 128 and we have chosen x = 3.15. Apart from 
the (expected) symmetry between the plotted curves, the figure shows that most of the monomers 
are located near the interface; the copolymer is thus adsorbed. Each of the two curves peaks near 
the interface, on the favorable side. The curves decay smoothly on the favorable side and sharply 
across the selective interface. There is evidence of frustration in the form of a secondary peak in 
monomer density on the unfavorable solvent side of the interface. This situation arises because 
the covalent bonding of the chain forces some monomers (whose neighbors are in majority of the 
other type) to be immersed in the "wrong" solvent. This situation is due to an energetic-entropic 
balance: configurations with monomers in the "wrong" solvent are energy-unfavorable but favor- 
able in view of entropy. 

Evidently, in Fig. 1 the asymmetry of each curve with respect to the interface is related to the 
value of X, a fact reminiscent of paramagnetism, where applying an external field increases (de- 
creases) the numbers of spins parallel (antiparallel) to it. Here we determine the magnitude of the 
asymmetry across the interface by evaluating 



H H/2 H H H/2 

M = ^|p^(z)-ps(2)| =^PaW- E pb{z)-J2pb{z). (1) 

2 = 2 = 2 = ff/2 + l 2 = H/2+l 2=0 

Figure 2 shows the relation between M and x for different copolymer lengths N. Note the very 
small dependence of M on N; the curves for different N almost coincide. As in paramagnetism, 
M is a linear function of x for small x, and reaches a constant (here unity due to normalization) 
for large x- 



We now turn to the question how the copolymer behaves around the interface. For this we compute 
Rgz, the z-componcnt of the radius of gyration of the copolymer. Here Rg^ — X]t=i (-^« ^ 'z)^ /N, 
where Zi is the z-component of the ith monomer's position and z = X^i^i Zi/N. We are interested 
in the dependence of Rgz on x ^nd hence display in Fig. 3 Rgz{x)/ Rgz{0) as a function of \fNx 
for different A^. Notice that (except for very large x) all data collapse into a single curve. This 
demonstrates that Rgz{x) scales with for small and moderately large x- Scaling fails when x 

is very large, because then the copolymers are practically squeezed on the interface, having the A- 
B covalent bonds at the interface. In the region in which the curves collapse we observe a constant 
regime for very small x followed by a power-law decay regime for moderate \/N x values. Setting 
Rgz{x)/Rgz{0) = fiVNx) with f{y) = 1 for j/ < and f{y) = y"" for y > yc, we find from Fig. 3 
that the crossover value t/c is roughly yc ~ 10. Furthermore, from a best fit to the data in the power- 
law regime we obtain numerically for the exponent a that a — 1.112 ±0.10 i.e., a/2 = 0.56 ±0.05. 
This result is in agreement with the scaling arguments of Refs. |l2j, which predict that a should 
equal 21/, where v is the usual Flory exponent, v « 0.588 for excluded-volumc chains. Scaling with 
ViV X is consistent with the "blob" -picture of a random copolymer at an interface |l^ . A blob is 
a chain segment containing g (g 3> 1) monomers, and thus has roughly g^^^ monomers of one type 
in excess; the total number of blobs is N/g and their number at the interface depends on \/N x [p^ . 



Another example of scaling is provided by the monomer density exactly at the interface z = 
[H + l)/2. In our (discrete) lattice model we thus compute 

Ps^\[pa{^) + Pb{^ + 1)] ■ (2) 

In Fig. 4 we display in a log-log plot ps as a function of x for different polymer lengths N. Again 
the data for different N coincide, possibly with the exception of the very small x regime; here, 



however, ps is small and thus the relative error is rather large. For very large Xj Ps gets to be 
independent of x, as the chains are then squeezed on the surface. In the moderate x regime ps 
scales with x as a power-law ps ^ i which is illustrated by the linearity of the data in the 
log-log plot of Fig. 4. From a best fit to the data in this regime we obtain numerically that the 
exponent f3 is 1.14 ± 0.06. This result can be compared to the expression fi = 2v in the "blob" 
picture |l^, i.e., taking the accepted v ~ 0.588 for excluded volume chains to (3 ~ 1.176. The 
derivation oi (3 = 2v according to Refs. |^ and jl^ starts from the scaling of physical variables 
with N^/'^x- Hence the total number of monomers at the interface obeys Ng — N^^"^ /(N^^'^x): 
where /(O) — 1. Noting that for large x is proportional to N requires that f{y) ^ y"\ with 
1 — z/+m/2 = 1. Hence m ^ 2v and Ng ~ Nx^'^ ■ Due to normalization one finds ps ~ Ng/N ~ x^". 

We conclude by showing that even the total monomer density scales with x'^^ ■ In Fig. 5(a) we plot 
{pa{z) + pb{z))/po as a function of (z- zo)x^'' with zq = {H+l)/2 and po = ^[pa(§) + Pb{§) + 
Pa{y + 1) + Ps(-f" + I)])- For a series of N and x values we find that all data collapse to a single 
curve when x lies in the power-law regime of Fig. 3. This implies the scaling relation 

{pa{z) + pb{z)) ^ pog{{z - zo)x^'") ■ (3) 

Eq. (3) does not hold when x lies outside the power-law regime. Fig. 5(b) displays the data of Fig. 
5(a) in a semi-logarithmic plot. The appearance of two almost straight lines rules out a Gaussian 
behavior and suggests to approximate the wings of g(y) by an exponential form. We find for the 
wings g{y) ~ exp {~j\y\), with 7 = 0.10 ± 0.02. 

In conclusion, we have studied numerically the adsorption of random copolymers at a selective 
interface. The main results are as follows. For small and for moderately large x the magnitude 
M of the asymmetry, as defined in Eq. (1), grows linearly with x and is independent of N. Fur- 
thermore, Rgz, the radius of gyration in the direction perpendicular to the interface scales as 
Rgz — f{-/Nx)- The monomer density ps at the interface, Eq. (3), scales as x^"; moreover, 
for X in the power-law regime the total monomer density obeys pa{z) + pb{z) = Pogi\z — zolx^") 
where g{y) is close to being exponential. 
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Figure Captions 

Figure 1: The probability (normalized density) p to find a monomer of A-type (diamonds) or 
B-type (crosses) at height z. Here the copolymer length is = 128 and x = 3.15. 

Figure 2: The dependence of the asymmetry parameter M, Eq.(l) on x for different polymer 

lengths. 

Figure 3: Rgz{x)/Rgz{0) as a function of the scaling variable VNx foi" different polymer lengths. 

Figure 4: The density ps exactly at the interface versus the interfacial selectivity strength x for 
different polymer lengths. 

Figure 5(a): The total monomer density plotted as a function of {z — zq)x^'' , see text for details. 
The parameter values are a: N = 32, x = 3.15; b: N = 64, x = 2.20; c: A = 64, x = 3.15; d: 
N = 128, X = 1-55; e: N ^ 128, x = 2.30; f:N= 128, x = 3.15. 
Figure 5(b): Same as (a) but plotted semi-logarithmically. 
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